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ABSTRACT 

Radial velocity measurements of stellar reflex motion have revealed many extrasolar planets, but 
gaps in the observations produce aliases, spurious frequencies that are frequently confused with the 
planets' orbital frequencies. In the case of Gl 581 d, the distinction between an alias and the true 
freq uency was the distinction betw een a frozen, dead planet and a planet possibly hospitable to 
life (jUdrv et al.l l2007t iMavor et al.l 12009). To improve the characterization of planetary systems, 
we describe how aliases originate and present a new approach for distinguishing between orbital 
frequencies and their aliases. Our approach harnesses features in the spectral window function to 
compare the amplitude and phase of predicted aliases with peaks present in the data. We apply it 
to confirm prior alias distinctions for the planets GJ 876 d and HD 75898 b. We find that the true 
periods of Gl 581 d and HD 73526 b/c remain ambiguous. We revise the periods of HD 156668 b and 
55 Cnc e, which were afflicted by daily aliases. For HD 156668 b, the correct period is 1.2699 days 
and minimum mass is (3.1 ± 0.4) M^. For 55 Cnc e, the correct period is 0.7365 days - the shortest 
of any known planet - and minimum mass is (8.3 ± 0.3) M^,. This revision produces a significantly 
improved 5-planet Keplerian fit for 55 Cnc, and a self-consistent dynamical fit describes the data just 
as well. As radial velocity techniques push to ever-smaller planets, often found in systems of multiple 
planets, distinguishing true periods from aliases will become increasingly important. 

Subject headings: planetary systems-planets and satellites: individual HD 156668 b-planets and satel- 
lites: individual 55 Cnc e-planets and satellites: individual GJ 876 d-methods: 
data analysis-techniques: radial velocities 
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1. INTRODUCTION 

In the past two decades, over 400 extrasolar planets 
have been discovered, including more than 300 detected 
by radial velocity measurements. The entire architecture 
of a planetary system is encoded in the wobbles of its host 
star. In frequency space, the star's radial velocity varia- 
tions are decomposed into the frequencies associated with 
each planet's gravitational interactions. One obstacle in 
correctly attributing these frequencies to planets are the 
spurious alias frequencies in the periodogram of the star's 
radial velocity measurements, caused by the discrete time 
sampling of the observations. Convolved with the orbital 
frequencies of alien worlds are Earth's own rotational and 
orbital frequencies, which dictate when the host star is 
visible at night, and - for many data sets - the synodic lu- 
nar frequency, which impacts the allocation of telescope 
time. 

Distinguishing aliases from physical frequencies is a 
common problem, yet making the correct distinction is 
cruc ial for characteriz ing extrasolar planets. For exam- 
ple, lUdrv et al] (|2007[ ) announced a super-Earth orbit- 
ing the M star Gl 581 with period 83 days, beyond the 
cold edge of the habitable zone. After more than dou- 
bling the number of observations, they determined that 
the planet's period was actually 67 days, well within 
the habit able zone, and th at the 83 day period was 
an alias (jMavor et al.l 120091 ) . The distinction between 
an alias and physical frequency was the distinction be- 
tween a frozen, dead planet and a planet possibly hos- 
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pitable to life. For reasons we will describe below, planets 
with periods of one to several months - in or near the 
habitable zone of M stars - will typically have aliases 
with periods within about 30 days of their own or- 
bital period. As more planets are discovered orbiting 
M stars, astronomers will be struggling to distinguish 
which of two close frequencies, one of which places the 
planet in the habitable zone, corresponds to a planet's 
orbital frequency. In general, planets with periods be- 
tween a few months and a few years often have confus- 
ing aliases caused by convolution with Earth's orbital 
period, while planets with periods near a d ay, such as 
the super-Earth G J 876 d (Rivera et al.ll2005l) , have con- 
fusing aliases caused by convolution with Earth's rota- 
tional pe riod. Automatic de- aliasing algorithms, such as 
CLEAN (jRoberts et al.lll987t) . have been applied to par- 
ticularly complicat ed radial velocity periodograms with 
some success (Oueloz et al.l l2009l ). yet, while they are 
good for cleaning up a periodogram, they should not be 
relied on for distinguishing between an alias and a phys- 
ical frequency. Aliases also pose a challenge for observ- 
ing variable stars and period-searching algorithms have 
been designed to not fall prey to them ( see for example 
iPlavchan et ^1120081: lReegenl[2007l Mm . 

Therefore, to enhance detection and characterization 
of planets, we have developed an approach to identify 
aliases by harnessing features of the "spectral window 
function," the Fourier transform of the observation times. 
Consider the star's motion as a signal that passes through 
a system, the time sampling window. Because of noise 
and loss of information, we can never perfectly recon- 
struct the signal. But we know everything there is to 
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know about the system: for a sinusoid of a given ampli- 
tude, frequency, and phase, peaks in the window func- 
tion cause ahases with calculable amplitudes and phases 
(|Deemingjll975l |1976[ ). The several time sampling fre- 
quencies - sidereal year, sidereal day, solar day, and syn- 
odic month - complicate the radial velocity periodogram 
yet allow us to break the degeneracy between alias and 
physical frequency that would exist for evenly-sampled 
data. 

In the following section, we describe the origin and 
characteristics of aliases, supply the details of our ap- 
proach for confirming that a particular frequency is 
not an alias, and clarify previous misconceptions about 
aliases. In the third section, we apply our approach to 
confirm periods for the planets GJ 876 d and HD 75898 b. 
We find that the orbital period for Gl 581 d and for the 
planets of HD 73526 cannot be definitively determined 
due to noise. We discover that the reported orbital pe- 
riod for HD 156668 b, 4.6455 days, is an alias of the true 
period, 1.2699 days. Finally, we analyze the 5-planet 
system 55 Cue. We find that the peri od of 2.817 days 
reported in the literature for planet e (jMcArthur et al.l 
l2004t iFischer et"al1 120081 ) is actually a daily ahas of its 
true period of 0.737 days. We conclude by summarizing 
the approach we have developed, considering the impli- 
cations of a new period for 55 Cnc e, and suggesting 
observational strategies for mitigating aliases. 

2. METHOD 

The existence of a planet orbiting a star is frequently 
inferred from a signature peak in the periodogram of ra- 
dial velocity measurements of the star. However, the 
periodogram often contains alias frequencies, the result 
of discrete sampling times, that, at first glance, can- 
not be distinguished from the true periodicities. Many 
astronomers have struggled to determine which peri- 
odogram peaks are physical frequencies and which are 
aliases, often resorting to methods that are unnecessarily 
computationally intensive, not definitive, refiect a misun- 
derstanding of aliases, or all of the above. In the first sub- 
section, we will describe the origin of aliases for evenly 
and unevenly sampled data. In the second subsection, 
we will explain the cause of the daily aliases, prominent 
for many Doppler datasets. In the third section, we will 
present a field guide for identifying aliases. In the fourth 
section, we will describe the method we have developed. 
In the fifth subsection, we will discuss the effects of or- 
bital eccentricity. In the sixth subsection, we will dis- 
cuss common misconceptions about aliases that lead to 
misidentification. 

2.1. The Origin of Aliases for Evenly and Unevenly 
Sampled Data 

Aliases are the result of discretely sampling a contin- 
uous signal. The resulting discretely-sampled signal is 
the product of the continuous signal and the sampling 
function, the latter being a "Dirac comb": a series of 
delta functions. The periodogram of the discretely- 
sampled signal is a convolution of periodogram of the 
continuous signal and the periodogram of the sampling 
function (the spectral window function). Consider first 
the simplified case of an infinite set of evenly-spaced 
data points, (?i[n], the result of sampling a continuous 
sine wave si{t) of frequency / at sampling frequency fs- 
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Here we follow iMcClellan et al.l (fl999h : 

si(t)=sin(2^/t), 
gi[n]=s{n/fs) = sin(27r/n//s), 

where n is an integer. However, under this sampling, 
the signal is indistinguishable from the sine wave S2 (t) of 
frequency {f + mfs): 
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where m is an integer. In the frequency domain, both gi 
and 32 will have peaks not only at /, but also at f + mfs. 
Moreover, neither has a periodogram distinguishable 
from a sampled sinusoid of frequency (— / -I- rnfg): 
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= sin(27r(-/ + mfs)n/ fs) = su\{2tt f n / f ^ + tt). 



That is to say, gi and g2 will also have peaks at -— / + 
mfs , although the phase of those peaks will be advanced 
by ^ cycle. For evenly sampled data, unless the only 
physically possible frequencies fall in a single Nyquist 
interval /i/2, the frequency cannot be unambiguously 
determined. 

Fig.[T]shows the spectral window function of an evenly 
sampled time series of /^ = 1 day~^. Peaks in the spec- 
tral window function occur at w./s, where m is an integer. 
The spectral window function is given by equation 8 in 
[Roberts et"all (|1987[ ): 
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where N is the number of data points, and t^ is the time 
of the rth data point. It is evident that when i/ = zLmfs, 

gT2^™/,t,, ^ gT2^™n ^ I ^^^ ^(j,) ^ I j^,g g^lgQ 

evident from this equation that when z/ = 0, W{i') = 1. 
Note that W{-i^) = W*{iy). 




Freq (1/day) 

Fig. 1. — Spectral window of data evenly sampled in time, with 
a sampling frequency /s = 1 day~^, and 300 samples. 



The top panel of Fig. [2] shows the periodogram^ of 
a sinusoid of period 1.94 days sampled every 1 day for 
300 days. For a sinusoidal signal, the resulting peri- 
odogram is a convolution of the spectral window func- 
tion W{i') with the peak corresponding to period 1.94 
days. The bottom panel shows the periodogram of a si- 
nusoid of period 2.06 days, an alias of 1.94 days, with 

^ For this and all other periodograms in this paper, at each fre- 
quency we (1) let the mean of the data float, and (2) weighted each 
dat a point with the i nverse of t he square of the reported error bar. 
See lCumminget al.1 l fl999f) and lZechmeister fc Kiirster) 120091) . 
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the same even sampling. The two periodograms are in- 
distinguishable. The aliases of the 1.94 day period occur 
at f — 1/1.94 + mfs- For fs~l and m — —1, the alias 
is i/1.94- 1 = 1/2.06. 

For a randomly selected frequency i^ each g^^'^*"*'' will 
add incoherently. However, if there are gaps in the data 
of a certain frequency ly, only certain phases occur and 
the complex exponentials will add in a partially coherent 
manner. The spectral window functions of stellar reflex 
motion measurements contain peaks at 1 sidereal year, 1 
sidereal day, 1 solar day, and sometimes 1 synodic month. 
These periodicities are caused by observations being lim- 
ited to only a particular portion of each of these periods. 
Observations are limited to a particular portion of the 
sidereal year and sidereal day because the star is only 
visible at night from the location of the telescope dur- 
ing particular parts of the sidereal year and day. At 
some telescopes, spectroscopic observations of the stars 
are relegated to "bright time," the portion of the synodic 
month when the moon is near full, because "dark time" 
is reserved for observing faint objects. In the next sec- 
tion, we will focus on the daily aliases due to both the 
solar day and the sidereal day. 

Uneven sampling also dict ates that the phase o f 
exp(27ri/sir) will span a width. lEver fc Bartholdl ()1999f ) 
demonstrate that for unevenly sampled data, there is ef- 
fectively no Nyquist frequency. Because gaps in the data 
and uneven spacing sample a non-zero width in phase, 
the height of peaks in the window function will never be 
exactly 1. For a noiseless data set, the physical frequency 
will almost always be a higher peak in the periodogram 
than any alias. (The only exception is if positive and neg- 
ative aliases add coherently.) For noisy data, the noise 
between two candidate peaks is correlated, but it may 
constructively interfere with the alias and destructively 
interfere with the true frequency, resulting in the alias 
peak being taller. Depending the phase of noise, it can 
also alter the phase of the true frequency and aliases 
through vector addition. 

2.2. Daily Aliases 

For most Doppler datasets, the largest peaks in the 
window function — corresponding to the largest aliases 
— are those at n day~^, where n is an integer. We refer 
to these peaks as the daily aliases, as they result from 
the sampling an Earth-bound observer is able to do at 
nighttime from a single site. 

Let us construct an example dataset, to illustrate their 
origin. Suppose the sampling is confined to when the 
Sun is down and the target star is up. In particular, sup- 
pose the samples are taken nearly daily, midway between 
when the star rises and the Sun rises, or midway between 
when the Sun sets and the star sets, depending on the 
time of the year. This sampling would lead to spacings 
between the solar day (24h Om Os) and the sidereal day 
(23h 56m 4s). Therefore, in our example dataset, let us 
take datapoints spaced by 23h 57m 30s, although due to 
telescope scheduling and weather, only a fraction of the 
nights (randomly chosen) are actually observed. Such a 
sequence is repeated in intervals of 365 days for 5 years, 
resulting in a total of 97 observation times. In Fig. |3] 
we illustrate this idealized dataset. It is constructed to 
obey the boundaries set by the Sun and the star, which 
are also plotted. The actual times from real datasets are 



compared, to show that this sampling, though idealized, 
reproduces the main daily and yearly structure of a real 
dataset. 

The window function for this idealized dataset is shown 
in Fig. m There are peaks at frequencies of n day^^ 



m yr 



In particular, there is a doublet at ly 



1.0000 day-i and j/ = 1.0027 day-^ with the latter peak 
being larger. 

How does this structure arise? We see from Fig. |3] that 
for 1/ = 1.0000 day~^, the idealized observations only 
sample the second half of phase. Therefore the window 
function as defined by equation ^ will have contribu- 
tions only from phases tt to 27r, so the complex exponen- 
tial will add up coherently to a large peak. This phase 
coherence explains the daily aliases not just at 1 day~^, 
but everywhere a peak occurs. For instance, consider 
at what times data are taken relative to the frequency 
of the sidereal day, v = 1.0027 day-^. In Fig. [3 this 
frequency is related to the diagonal line labeled "star 
rises." The idealized dataset consists only of observa- 
tion times between 0.1 days and 0.4 days after the star 
rises (above that diagonal line). Therefore the observa- 
tions cover only 30% of the phase of the sidereal sampling 
frequency, which again results in a large peak in the win- 
dow function. Here, even a smaller fraction of the total 
phase is covered, so the sampling results in even more 
coherent summation of complex exponentials, which is 
why the window function peak &t v = 1.0027 day~^ is 
larger than a± v ~ 1.0000 day^^ (Fig. |4|). Another way 
to see this is to note that the line formed by the idealized 
data in Fig. [31 panel a, has a slope more closely match- 
ing the sidereal day (the diagonal lines related to the 
star) than the solar day (the horizontal lines related to 
the Sun). Finally, we note that no peaks in the window 
function appear between the solar and sidereal frequen- 
cies because folding the data at those frequencies samples 
phases throughout to 27r. 

Having understood the origin of the daily aliases in the 
window function, including doublets, we are prepared to 
recognize and correctly interpret such structure when it 
results in periodograms. 

To that end, we used this idealized dataset to sample 
a sinusoid of period 1.94 days or 2.06 days, and in Fig. [5] 
show their periodograms. In this exa mple, we have taken 
the two periods close to those which ' Rivera et al.l ()2005f ) 
needed to decide between for GJ 876 d. Here, then, we 
have identified a simple way to decide between them: 
the slightly taller peak is expected to be the true one 
(because there is no noise), and the alias will consist of 
a doublet wi th spacing 0.0027 day~^. We analyze the 
iRivera et al.l ([2005) dataset in subsection 13. II 

2.3. A Field Guide to Aliases 

An alias is a convolution in frequency space of a physi- 
cal frequency with the window function. Fig.lHland Fig. [7] 
display some examples of yearly and daily aliases respec- 
tively and the window function features that cause them. 
We have chosen especially clean examples; ambiguous 
cases will be addressed throughout the next section. 

2.4. Details of Our Method 

We recommend the following treatment for a radial 
velocity dataset or residuals of an established fit (we 
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Fig. 2. — Periodogram of sinusoids sampled evenly in time. Top: Period 1.94 day. Bottom: Period 2.06 day. They are indistinguishable. 



will refer to both these categories as "data") that ap- 
pear to exhibit periodic variation. As we emphasized 
above, the phases of peaks are helpful for determining 
what is the true frequency. For example, consider a set 
of data with peaks in the spectral window function at 
1 year (0.0027 day^^), 1 solar day (1 day~^), and 1 
sidereal day (1.0027 day""'^). Consider a true frequency 
/i > 1.0027, which will have aliases at at /2 = /i — 1 
and /a = /i — 1.0027. We may wonder if the peak at /2 
is the true frequency, with an alias at fs — 0.0027 = /a. 
However, because of the phases of the peaks in the win- 
dow function, the phase of the peak f^ is different than 
the phase we would expect if it were an alias of /2 . Be- 
cause the phase of a peak can be key in determining 
the true frequency, we strongly recommend plotting the 
phase of selected peaks. We use a symbol we call a "dial" 
(e.g.. Fig. [9]) where the phase angle is the counterclock- 
wise angular position from the x-axis. The phase angle 
is tan~^(Imaginary(M^(i/))/Real(VF(:^))) for the window 
function peaks and likewise tan~^(C(/)/i?(/)) for the 
periodogram peaks, where B and C are the real and 
imaginary coefficients of the periodogram for frequency 

/• 
Our method is composed of the following steps: 

1. Plot the spectral window function (eq. [T]), attach- 
ing dials to any large peaks. Peaks will most likely 
occur at or near f=0, 1/yr, l/(solar day), 1/(1 
sidereal day), and, if the observations were taken 
during a particular part of the lunar cycle, f — 
1/month, 1/month ± 1/yr. Spectral window func- 
tions of artificial data sets are plotted in Fig. [T] and 
12 and real data sets in Fig. [HI UHl HH [HI US [HI [H 
and [211 

2. Plot the periodogram. 

3. Consider first the possibility that the largest peak 
is the true frequency; measure its frequency, phase, 
and amplitude. Attach dials to peaks we would 
expect are aliases, according to the peaks in the 
window function. If the peak in the radial velocity 
periodogram occurs at / and peaks in the window 
function occur at fs, we expect aliases at /zb/^. (If 
fs > /> we will still see a peak at f — fs- Flipping 
it across frequency gives the phase the opposite 
sign: a complex conjugation.) Generate a sinusoid 
with the same frequency, phase, and amplitude as 
the peak in radial velocity periodogram and plot its 



periodogram, attaching dials in the same location. 
Compare the amplitude and phases of peaks. Are 
the major aliases for / present in the data with the 
predicted phase and amplitude? 

4. Now consider that the largest alias(es) of what we 
considered the true frequency might actually be the 
true frequency. Repeat step 3. 

5. If the periodogram of the data is well-matched by 
the periodogram of one and only one candidate 
sinusoid, t hen the true frequency has been deter- 
mined. As iLombl ()1976l ) said, "If there is a satis- 
factory match between an observed spectrum and 
a noise-free spectrum of period P, then P is the 
true period." However, if several candidate sinu- 
soids match peaks equally well or poorly, then the 
data are not sufficient to distinguish the true pe- 
riod. 

2.5. Treating the Orbital Eccentricity 

Many extra-solar planets have elliptical orbits. The 
signal of the eccentricity is contained in harmonics of 
the orbital frequency; the first harmonic has an am- 
plitude eK, where e is the eccentricity and K is the 
amplitu de of the sine wave at the p lanet's orbital fre- 
quency (jAnglada-Escude et al.l [20101 ). Thus for moder- 
ate eccentricities, the same analysis can be applied to 
the first harmonic of the orbital frequency. Except in 
rare unfortunate cases (such as HD 73526, treated be- 
low in section 13. 3p . the period and its aliases will be 
well-separated in frequency space from the eccentricity 
harmonic and its aliases. In section [321 we distinguish a 
pea k in the peri odogram of HD 75898 b as an alias that 
the iRobinson et al. (2007) proposed could be an alias, 
eccentricity harmonic, or additional planet. 

For certain datasets, orbital eccentricity may help dis- 
tinguish between a true orbital period and an alias. Con- 
sider a planet with moderate eccentricity e whose host 
star is observed with near evenly-spaced sampling as fs- 
Even if the noise is low relative to K, it may be difficult 
to distinguish between the true orbital frequency / and 
an alias / -I- fs. However, since the planet's orbit is ec- 
centric, we will also observe a peak of amplitude eK at 
2/ but no such peak at 2{f + fs). 

In summary, orbital eccentricity contributes to the pe- 
riodogram in a well-defined way and, except in rare un- 
fortunate cases that can be easily identified, will not con- 
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Fig. 3. — Times of observation of an idealized dataset and two 
real datasets, folded to illustrate the origin of daily aliases. The 
axes show, quantitatively, the time of the year and the time of the 
day. The solid lines are labeled and correspond to the time each 
day that either the Sun rises or sets (at a constant time-of-day 
in this idealized example) or the star rises or sets (which varies 
according to the time of the year). The dashed lines are when the 
star reaches 54 degrees from the zenith, within which a favorable 
observation can be made. The idealized dataset is described in the 
text. The HARPS data for Gl 581 are from Mayor et al. (200^), 
and we took t = JD - 2, 452, 970.9 2 for convenience. T he Keck 
and Lick data for 55 Cnc are from IFischer et al.l l|2008l ') , and we 
took t = JD-2, 447, 370.15. 

fuse the distinction between the true orbital period and 
an alias. 

2.6. Common Misconceptions 

Many problems with aliases are the result of unwar- 
ranted assumptions. We describe some common miscon- 
ceptions about aliases and how they cause confusion. 

1. Assuming that the largest peak in the periodogram 
is the physical frequency. In fact, noise may add 
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Fig. 4. — Spectral window function of data with gaps. The 
sampling is from the "idealized dataset" of panel (a) in Fig. [3] 

coherently to an alias or incoherently to the physi- 
cal frequency, causing the alias to ap pear larger. 
This is what happened for Gl 581 (jUdrv et all 
120071 : iMavor et alJl2009l ). In muhi-planet systems, 
aliases from several planets cou ld add to ma ke the 
highest peak a spurious signal (|Fosteilll995[ ). 

2. Assuming that the frequency that yields the best 
Keplerian or Newtonian planet fit is the true fre- 
quency. As we saw for Gl 581 d, this is not always 
the case, due to noise. 

3. Assuming that aliases occur at frequencies only oc- 
cur near peaks in the spectral window function. 
We have seen authors plot the spectral window 
function below the periodogram of the data and as- 
sume that if a frequency in the data periodogram 
is not near a peak in the spectral window func- 
tion that it is not an alias. In fact, aliases occur 
at 1/ ± /s|, where fs is a feature in the spectral 
window function. Depending on the relative val- 
ues of / and fs, the alias might be anywhere in 
the periodogram. However, periodograms will con- 
tain peaks at the sampling frequencies if there are 
systemics linked with the observing pattern or if 
the peaks are aliases of a very low frequency sig- 
nal. We emphasize the difference between these 
two types of signals: the former is spurious and 
and the latter has an extra-solar origin but wrong 
frequency. We also emphasize the importance of 
employing the spectral window function to identify 
all major aliases, not just aliases or other spurious 
frequencies that occur at the sampling frequencies. 

4. Assuming that any frequency above 1 is an alias. 
As we mentioned above, there is effectively no 
Nyquist frequency for unevenly sampled data. 
Many authors cut off their periodograms at 1 
day~^, potentially missing out on or misinterpret- 
ing planets with orbital periods less than a day. We 
know such planets exist because they have been de- 
tected by transits. Moreover, because long period 
planets will have aliases near 1 day, a planet with 
orbital periods near 1 day is vu lnerable to being 
discarded as an alias (|Kanell2007l ). 
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— Top: Periodogram of sinusoid of period 1.94 day (frequency 0.515 day ^) with the ideahzed time sampling from Fig.|3] Bottom: 
.06 day (frequency 0.485 day~^). With this time sampling, the periods are distinguishable by the imprinting of the window function 
from Fig. [4] at / ± /s where / is the frequency of the sinusoid and fs of the window function feature. 
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Fig. 6. — Illustrative examples of yearly aliases taken from GJ 876 (top), 55 Cnc (middle), and HD 156668 (bottom), 
function is plotted on the left and the periodogram of the data near the candidate frequency on the right. The arrow in 
indicates the peak in the window function near 1/yr and the arrows in the right plots indicate the predicted location of the 
caused by this window function feature. 
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Assuming that aliases are so pernicious that one 
can never identify the correct period and should 
thus just pick the most sensible period. In fact, 
our method allows one to determine either a cor- 
rect period or that noise prevents the identification 
of the correct period. In the latter case, further 
observations should allow for a definitive determi- 
nation in the future. It is unwise to judge a priori 
which period is the most "sensible" period; as men- 
tioned above, planets have been found with periods 
less than a day. 

Assuming that if an alias frequency is used in a Ke- 
plerian or Newtonian planet fit, a peak correspond- 
ing to the true frequency will appear in residuals. 
This would only happen if the peak at the alias fre- 
quency is much smaller than the peak at the true 
frequency, relative to the noise. 



7. Assuming that if a frequency is an alias, it will 
appear in a periodogram of the data scrambled. 
Aliases are not caused solely by the spacing obser- 
vations; they are convolution of the spectral win- 
dow function with the periodogram of the data. 
Scrambling the data removes the true frequency 
and thus also removes the alias. 

8. Assuming that if you "fold" (i.e., phase) the data 
with a candidate period, a coherent pattern will 
emerge only if the candidate period is the physical 
period. In fact, a large alias, by its very definition, 
will also produce a coherent pattern. 

Another method we have seen applied to distinguish 
between two frequencies, one of which is an alias, is to 
generate thousands of mock data sets for each frequency 
by combining a sinusoid with simulated noise and then 
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Fig. 7. — Illustrative examples of yearly aliases taken from HD 75898 (top), GJ 876 (second row), 55 Cnc Fischer et al. data set (third 
row), and 55 Cnc combined data set (bottom row). The window function near a major feature is plotted in the left column and sections of 
the periodogram of the data in the middle and right columns. Arrows in the left column indicate the peaks in the window function near 
sidereal and solar days and the arrows in the middle and right plots indicate the predicted locations of the corresponding aliases. Note 
that each peak in the window function results in two features in the data periodogram. 



determine how often the ahas is mistaken for the true 
frequency. This method indeed reveals the probability 
that the period is falsely determined, but a proper un- 
derstanding of the window function leads to a less com- 
putationally intensive method, which we have advocated. 
We reemphasize the peaks in the spectral window func- 
tions combined with the true frequencies are what cause 
aliases. Even if a peak in the periodogram is linked to an- 
other peak by close to an integer frequency, if that integer 
frequency is not a peak in the spectral window function, 
then the peaks are not aliases of one another and might 
represent two distinct planets. Rather than simply not- 
ing the possibility that an integer frequency might link 
the peaks, the window function reveals it quantitatively. 

3. APPLICATION TO EXTRASOLAR PLANETARY 
SYSTEMS 

In the following section, we investigate instances of 
aliases and ambiguous periods in the literature. 

3.1. GJ 876 d 

In this section, we apply the approach described above 
to planetary system GJ 876. Extensive radial velocity 
observations spanning almost eight years have revealed 
three planets orbiting this M-star. A Jupiter -mass planet 
b was discovered in 1998 (jMarcv et al.lll998l) . and an in- 
terior Jupiter-mass planet c in a 2:1 resonance w ith b 
was discovered three years later (jMarcv et al.ll200lD. Af- 
ter several years of continued observations, iRivera et al.l 



(|2005f ) discovered an additional 7.5 earth mass planet d 
with an orbital period of 1.9 4 days. This d i scover y was 
independently confirmed by iCorreia et al.l (|201C1[) with 
new HARPS data. The periodogram of the residuals 
to the nominal two-planet, i = 90°, coplanar fit ex- 
hibits strong power at frequency 0.52 day"""^ but also at 
/ = 0.49 day-i (P = 2.05 day) a nd / = 1.52 day"^ 
[P = 0.66 day) (Fig. El top panel). iRivera et al.1 ()2005f ) 
performed a series of tests and argued based on the re- 
sults that the peak at 2.05 days is an alias of the true 
period at 1.94 days. Our method is able to definitively 
confirm the results of their tests, that the physical period 
is indeed 1.94 days. 

The spectral window function and the periodogram of 
GJ 876 (actually, of the residuals from a dynamical fit to 
planets b and c) are shown in Fig. |8]and|9l respectively. 
Major peaks in the window function occur at 1 sidereal 
year, 1 sidereal day, and 1 solar day. The very same 
features are seen in the example periodogram described 
in section 12.21 The main peak is tallest "^ . The alias has 
a doublet structure. 



* We point this out for identification purposes but in a given 
data set, because of noise, the true frequency will not necessarily 
be taller than the alias. 



Dawson & Fabrycky 



Compared to our example idealized data set of Fig. [3l 
m and [5l the yearly aliases are more pronounced in the 
data, because the observing season is shorter than in 
our idealized dataset. This causes peaks on either side 
of the true peak (spaced by 1 yr^^ ~ 0.0027 days~^) 
which are symmetric in height. Thus we confirm the se- 
lection of P = 1.94 days as the correct period of GJ 876 d 
(jRivera et al.ll2005[ ) , and thus we demonstrate that a sig- 
nal beyond the traditional Nyquist frequency can be ro- 
bustly detected with unevenly sampled data. 
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Fig. 8. — Spectral wind ow function of RV Measurements of 
GJ 876 URivera et al.l [20051) . Major features of the spectral win- 
dow function are colored: red (at day~^), green (yearly fea- 
ture), fuscliia (daily features), blue (two day~^), and brown (three 
day~^). The corresponding aliases these features cause for several 
candidate frequencies are indicated by these colors in Fig. O 



3.2. HD 75898 b 

iRobinson et al.l (|2007[ ) discovered a Jupiter-mass 
planet orbiting HD 75898 b. They noticed two peaks in 
the periodogram, a large one near 400 days and a smaller 
one near 200 days. They presented three possibilities for 
the peak near 200 days: an alias of the 400 day period, 
an eccentricity harmonic (which we would indeed expect 
to appear near P/2 = 200 days), or a second planet. Ap- 
plying our method, we confirm that the true period is 400 
days, not 200 days; and the peak at 200 days is indeed an 
alias, not an eccentricity harmonic or second planet. The 
spectral window function is plotted in Fig. [TU] the peak 
that occurs at 1 yr~-^ is the cause of the 200 day alias. 
In Fig. [m the periodogram shows that a 400 day period 
(row 2) produces exactly the aliases we expect, including 
the alias at 200 days. Although an eccentricity harmonic 
would fall at the same place as this alias, for this system 
we can rule out a significant eccentricity harmonic be- 
cause the peak has the exact phase and amplitude that 
result from it being an alias of the 400 day planet; any 
significant eccentricity harmonic would change the phase 
and/or amplitude of this peak. These plots also confirm 
that the true period is 400 days, not 200 days (row 3). 

3.3. HD 73526 

iTinnev et al.l (|2003[ ) reported a planet orbiting the 
G-type star HD 73526 with orbital period 190.5 days. 



A later Bayesian analysis by iGregorvl ()2005l) revealed 
three possible periods for the planet: 190.4 days and 
(its yearly aliases) 127.88 days and 376.2 days. IGregorvl 
(2003) concluded that the periods 127.88 days and 376.2 
days were mo re probable. After follow-up observations, 

iTinnev et al.l (p006) reported the system actually con- 
tained two planets, with orbital periods 187.5 and 376.9 
days, locked in a 2:1 resonance. The Keplerian fit using 
these two periods is an excellent match to the data, with 
ixlY^'^ — 1-09, but the dynamical fit for the system is 
substantially worse, with (xtY^'^ = 1-57. This implies 
that, though these periodicities may be strongly present 
in the system, the physical model of two planets orbiting 
with this period may need modification. Further com- 
plicating the interpretation of the system's periodicities 
is the degeneracy between the outer planet's eccentricity 
and the inner planet's mass - or even its very existence 
(jAnglada-Escude et al.ll2010t) . The window function for 
this system and a periodogram is plotted in Fig. |T2l The 
Keplerian fit has eccentricities of 0.4 for both planets, 

g essentially tuning the phase of the power at 187.5 days 
(the first eccentricity harmonic of 376.9 days) to account 
for both a possible planet there and aliasing from 376.9 
days; and introducing power at 93.8 days (the first eccen- 
tricity harmonic of 187.5 days and also a yearly alias of 
127 days). However th e eccentricities for the dynamical 
fit (jTinnev et al.l l2006[ ) are substantially lower, imply- 
ing that high eccentricities would cause dynamical inter- 
actions inconsistent with the data. It is possible that 
the periods 127.88 days and 376.2 days are incorrect but 
that by introducing a large eccentricity harmonic, the 
combination of orbital periods, eccentricity harmonics, 
and aliases match the periodicities of the data, which 
may be the result of different physical orbital frequen- 
cies. This system is complicated because of the degener- 
acy in frequency between resonant planets, eccentricity, 
and aliases. We recommend further observations and 
modeling of this system to confirm the orbital periods. 

3.4. Gl 581 d 

HARPS measurements have revealed four plan- 
ets orbiting the M dwarf Gl 581: a -- 2M© 
planet e (^ avor e t al.l l2009t ). Neptune- mass planet b 
(Bonfils ct al. 2005), and super-Earth planets c and d 
(Udrv ct al. 2007). Planet d was originally reported to 
have a period of 83 days, beyond the cold edge of the 
habitable zone. After further observations, the HARPS 
team announced that the true period of planet d is 67 
days, placing it within the habitable zone, and that the 
original 83 day period was a one year alias of the true 
67 day peri od. In Fig. [131 we plot the spectral window 
function of iMavor et al.l ()2009f) 's new data set. Promi- 
nent peaks are evident at 1 year, 1 sidereal day, and 1 
solar day. A periodogram of the data, with planets b 
and c subtracted (subtracting planet e made no signif- 
icant difference) and sinusoids of several candidate fre- 
quencies are plotted in Fig. jTH In the original data set, 
the highest peak in the periodogram was at 0.0122 day~^ 
(corresponding to a period of 83 days) . In the new data 
set, the highest peak is at 0.9877 day^^. The second 
highest pe ak is at 0. 0150 day^^ (67 days), the period re- 
ported bv IMavor et al.. (2009,). The 0.0122 day^^ peak 
and 0.0150 day~^ are linked by a feature in the window 
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Fig. 9. — Periodograms of GJ 876. The top row is the periodogram of the data. The second and third rows show the periodograms of 
sinusoids sampled at the times of the real data sets as solid lines; they also repeat the periodogram of the data as a gray background, for 
comparison. Dials above the peaks show the phase at each peak. Colors correspond to the feature in the window function that creates the 
particular alias (see Fig. Is}, with red being the candidate frequency, the green sidebands yearly aliases, and the fuschia, blue, and brown 
peaks daily, two day~^, and three day~^ aliases respectively. The second row is the periodogram of an injected sinusoid of period 1.94 
days (frequency 0.516 day~^). The third row is the periodogram of an injected sinusoid of period 2.05 days (frequency 0.487 day~^). The 
sinusoid of period 1.94 days matches the heights and phases of the peaks much better, both for the yearly aliases on either side of the 
main peak in column 2 and the daily aliases in the other columns. The two candidate frequencies have different types of aliases at different 
locations, allowing us to break the degeneracy. 
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Fig. 10. — Spectral window function of RV Measurements of 
HD 75898. These features, convolved with a planet's orbital fre- 
quency, cause the aliases evident in the periodogram in Fig. 1111 



function at 1 sidereal year. Yet neither produces an alias 
that corresponds to the other frequency with a phase and 
amplitude that match the data (first column of rows 2 
and 3). The highest peak, 0.9877 day~^, is linked to the 
peaks at 0.0122 day~^ and 0.0150 days~^ by the win- 
dow function peaks at 1 solar day and 1 sidereal day 
respectively; it better matches the phase and amplitude 
at these frequencies (row 4, column 1). This dataset has 
sampling which is too regular (Fig. [5]d), which resulted 
in pernicious daily aliases. However, there are discrep- 
ancies between the phase and amplitude of the aliases 
predicted by all three candidate frequencies. For exam- 
ple, at 1.99 days (column 4), the larger alias predicted 
for 0.9877 day~^ (linked by the large 1 sidereal day alias) 
is consistent in amplitude with the data while the other 
frequencies (linked by the smaller window function fea- 
ture at 2 days) predict aliases that are too small; how- 
ever, the phase for the 0.9877 day~^ alias is a bit off. 
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Fig. 11. — Periodograms of HD 75898. Dials above the peaks 
denote their phase. Row 1 shows the data. The other rows show 
sinusoids sampled at the times of the real data sets (solid line and 
dial), as well as the data again for reference (in gray). In Row 2 
the solid line shows, for these time samplings, the periodogram of 
a sinusoid of frequency 0.00236 day~^. For Row 3, it is for 0.00519 
day~^. We confirm that the peak at 0.00519 day~^ is an alias, not 
a second planet or eccentricity harmonic. In Rows 2 and 3, each 
peak results from the convolution of the sinusoidal frequency with 
the features in the spectral window function in Fig. 1101 



Although none of the frequencies is fully consistent, we 
shghtly prefer 0.9877 day-^, followed by 0.0150 day"! 
and 0.01 22 day~^. Howe ver, using the previous data 
set from lUdrv et al.l (|2007[ ). we favor (in order): 0.0122 
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Fig. 12. — Top Panel: Spectral window function of RV mea- 
surements of HD 73526. Bottom Panel: Periodogram of RV mea- 
surements of HD 73526. The solid arrows indicate the locations 
of a peak's yearly aliases and the dashed line the location of the 
eccentricity harmonic. 

day-\ 0.0150 day-\ and 0.9877 day^ We also fit a 
four-planet Keple rian model to both datasets. In the 
lUdrv et all (|2007D dataset, a frequency of 0.0122 day"^ 
for pla net d gave the best fit, while in the iMavor et alJ 
(|2009t) dataset, a period of 0.9877 day-^ gave the best fit. 
However, a model with orbital frequency 0.0122 day~^ 
where e^ is allowed to float gives a significantly better 
fit than one with orbital frequency 0.9877 day~^ where 
6(1 is fixed at zero (which would likely be attained by 
tidal dissipation). Because the period of planet d remains 
ambiguous, we recommend that future observations take 
place with the star at a greater air mass - instead of only 
when the star is crossing the meridian - in order to reduce 
the amplitude of the aliases and allow us to definitively 
distinguish between these three candidate periods. 

3.5. HD 156668 b 



iHoward et al.l (|2010( ) reported a 4 M^ planet orbiting 
HD 156668 b with period 4.6455 days (a frequency of 
0.2153 day~^). However, they considered that the cor- 
rect period might be 1.2699 days (a frequency of 0.7875 
day~^), and our analysis confirms that as the correct pe- 
riod, as follows. The window function for this system is 
plotted in Fig.fTSland periodograms of the data and sinu- 
soids at two candidate frequencies in Fig. [161 Note that 
large peaks in the window function occur at 1 sidereal 
and 1 synodic day while smaller peaks occur near 2 days 
(Fig. [T5|). For a true frequency of 0.2153 day~^ (second 
row), we would expect two pairs of large peaks due to 
sidereal and solar aliases (second row, second and third 
column) and a smaller pair of peaks for the ^2 day^^ 
aliases (second row, fourth column). On the other hand, 
for a true frequency of 0.7875 day^^ (third row), we 
would expect two pairs of large peaks due to sidereal and 
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Fig. 13. — Spectral window function of Gl 581. These features, 
convolved with a planet's orbit al fr equency, cause the aliases evi- 
dent in the periodogram in Fig. 1141 

solar aliases (third row, first and fourth column) and a 
smaller pair of peaks for the ~ 2 day~^ aliases (third row, 
third column) . The phase and amplitude of these aliases 
predicted for 0.7875 day~^ (row 3) are thus more consis- 
tent with the data (row 1). Therefore we conclude that 
the planet's true period is 1.2699 d ays and that the peak 
at period 4.6455 days identified bv lHoward et al.l ()2010f ) 
is an alias. The Keplerian orbital elements are reported 
in Table [H along with the predicted transit window. The 
eccentricity was he ld to zero, as expecte d from tidal dis- 
sipation, following iHoward et al.l (|2010l ). iHoward et al.l 
(2010) "filtered" the data by simultaneously fitting a 
two-planet model and a linear trend. They state that 
the "second planet" is a form of high-pass filter, not nec- 
essarily an actual planet. We do not fit a linear trend or 
additional planets in our reported fit and do not subtract 
them out in Fig. 1161 However, we have confirmed that 
our results hold if we do. 

3.6. 55 Cnc 



With five discovered planets (jFischer et al.l2008[ l , more 
than any other extrasolar planetary system, 55 Cnc is a 
rich en vironment for study. The first planet was discov- 
ered bv lButler et al.l (|1997| ): this planet b has an orbital 
period of 14.65 days. Five mo re years of observa tions re- 
vealed two additional planets ()Marcv et al.ll2002l ): planet 
c, with orbital period 44 days, and planet d, with orbital 
period 500 days. Measurements f rom the Hobby- Eberly 
Telescope (jMcArthur et al.l l2004l ) (HET) revealed, on 
th eir own and c ombin ed with the Lick measurements 
bv iMarcv et al.l (l2002t ) and ELODIE measurements by 
iNaef et all (|2004l ). the presence of planet e, with a re- 
ported orbit al period of 2.8 days. In 2005, in a poster 
presentation (|Wisdomll2005l) and an informally circulated 
paper ^, Wisdom (hereby referred to as W05) reanalyzed 
the combined HET, Lick, and ELODIE measurements, 
found evidence for a 260 day period planet, and ques- 
tioned whether the repor ted 2.8 day signal m ight be an 
alias of planet c. Finallv. [Fischer et al.l ( 2008| ) confirmed 
the 2.8 day planet e and reported a 260 day planet f based 
on a decade of Lick and Keck measurements. They also 
noted a peak at 460 days and considered whether this 

^ Available electronically at 

[http://groups.csail.mit. edu/mac/users/wisdom/planet.ps^ 
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Fig. 14. — Periodograms of Gl 581 for planet d (planets b and c subtracted have been removed from the data set and planet e has been 
ignored; we obtain consistent results if we also remove planet e). Dials above the peaks denote their phase. Row 1 shows the data. The 
other rows show sinusoids sampled at the times of the real data sets (solid line and dial), as well as the data again for reference (in gray). 
In Row 2 the solid line shows, for these time samplings, the periodogram of a sinusoid of frequency 0.0122 day^^. For Row 3, it is for 
0.0150 day~^. For Row 4, 0.9877 day^^. In Rows 2-4, each peak results from the convolution of the sinusoidal frequency with the features 
in the spectral window function in Fig. 1131 Note that the phases and amplitudes of 0.0122 day~^ and 0.0150 day~^ are not consistent 
with the aliases we would expect. The period remains ambiguous, but we favor 0.9877 day~^ based on this data set. 



TABLE 1 
New parameters for HD 156668, e^ = 0.' 



K 

ms~ 



M sin i 



p 

days 



AU 



deg 



A 
deg 



V 



2.2(3) 3.1(4) 



1.26984(7) 0.0211(2) 0.000(0) 0.(0) 



136.(19) 



-0.4(2) 



Note. — Data arc the Keck data presented by lHoward et al.j <2Q1(]| ). T^ ^^-^ is set to the first data point (JD 2453478.97768). These parameters prediet a transit epoeh 
of TtrfJD] = 2453478.82(7) + E X 1.26984(7). 

^ The foUowing gravitational constants were used: GMq = 0.0002959122082856, ratio of the sun to Earth = 332945.51. The mass of the star was assumed to be 0.77 solar 

^iven in parentlreses, referring to tire final digit(s). 
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Fig. 15. — Spectral window function of RV measurements of 
HD 156668. Major features of the spectral window function are 
colored: red (at day~^), green (yearly feature), fuschia (daily 
features), and blue (two day~^). The corresponding aliases these 
features cause for several candidate frequencies are indicated by 
these colors in ll6l 



peak was an alias of the 260 day planet. 

Because the literature has considered whether they 
might be aliases and because their periods are in the 
range where aliases can be the most confusing, planet e 
and planet f warrant additional consideration. We con- 



firmed by our analysis that the period of f is correct. In 
the following subsection, we apply our method to planet 
e and find that the 2.8 day period is actually an alias, 
not of planet c but of a true period of 0.74 days: planet 
e still exists but its period is actually 0.7 days, not 2.8 
days. 

3.6.1. A New Period for 55 Cnc e 

First, let us look at the discovery data for 55 Cnc e. 
We plot the window f unction for the data collected by 
iMcArthur et all ()2004D using HET in Fig. HT] The data 
spans only 190 days and therefore contains no yearly 
gaps. Therefore, no peak in the window function oc- 
curs at yr~^, and there is no splitting of the daily alias 
into solar and sidereal days. We also note that this daily 
alias has quite a strong value of ~ 0.8. The consequence 
of that can be seen in Fig. [HI the periodogram using 
only the HET data. The top panels are the periodogram 
of the data themselves. The peaks at 0.356 day~^ and 
1.358 day^^ are of similar size. In the middle panels, we 
sample a noiseless sinusoid with a period, amplitude, and 
phase matching that of the peak at 0.356 day~^. An alias 
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Fig. 16. — Pcriodograms of HD 156668. Row 1 shows the data. 
The other rows show sinusoids sampled at the times of the real data 
sets (solid line and dial), as well as the data again for reference (in 
gray). Colors correspond to the feature in the window function 
that creates the particular alias (see Fig. I15I I. with red being the 
candidate frequency, the green sidebands yearly aliases, and the 
fuschia and blue peaks daily and two day~^ aliases respectively. 
In Row 2 the sohd hne shows, for these time samplings, the peri- 
odogram of a sinusoid of frequency 0.215 day~^. For Row 3, it is 
for 0.787 day~^, our favored value. The two candidate frequencies 
have different types of aliases at different locations, allowing us to 
break the degeneracy. 



results at 1.3 58 day~^ at approxima tely the right height 
and phase, so lMcArthur et all ()2004[ ) may have dismissed 
the latter as an alias, although they did not mention it 
explicitly. However, reversing the argument, if we had a 
noiseless sinusoid with the period, amplitude, and phase 
of the peak at 1.358 day~^ (bottom panels), then its alias 
nearly matches the peak a t 0.356 day~^, wi t hin th e noise. 
This is to say, the data of iMcArthur et all (|2004[ ) cannot 
distinguish between the two possible periods. 

W05 presented two arguments for why the 2.8 day sig- 
nal might be an alias. First he noticed that the 2.8 day 
period is linked to the 44 day period of planet c by a 
period of 3 days (2^ w 5 + s)' ^^"^^ noted that there 
is no reason we would expect an alias to be caused by 
a 3 day period. In Fig. [171 US and [21] we demonstrate 
that there is no peak in the spectral window function at 
i day^^ for any of the data sets. Therefore, the 2.8 day 
signal cannot be an alias of the 44 day signal. 

Second, W05 noticed that in the HET data, one peak 
occurs at 2.808 day, while in the combined data set a pair 
of peaks occurs at 2.7957 days and 2.8175 day, a split- 
ting of 1 year. In fact, this is just the doublet structure 
described in section 12.21 The combined set spans multi- 
ple years, which creates the yr~^ spacing in the doublet 
structure of the daily alias, as shown in Fig. [THl There- 
fore we would actually expect to see this doublet struc- 
ture in the combined data set but only a single peak at 
the daily aliases in the HET data set. 

So Wisdom was right to suspect that the 2.8 day signal 
is an alias. It is not an alias of the 44 day planet c but 
of a planet with true period 0.7 days; the alias is a daily 
alias (1/2.8 days = 1/0.74 days - 1/days). 

With the combined data set, and with new data that 
has come out with higher precision from Lick and Keck 
([Fischer et al.l i2008.) , we can confirm with high confi- 
dence that the 0.74 day period is the correct one. The 



window functions of these datasets are shown in Fig. [191 
and [211 In Fig. [201 and HH we show the resulting pcri- 
odograms, after subtracting the signal of planets b, c, and 
d with a best-fitting Keplerian model. In both datasets, 
the true peak at 1.358 day~^ is very much higher and 
the other peaks at various frequencies are fully consistent 
with being an alias of it. For instance, in both datasets, 
doublet structure at the reported frequency shows that 
it is actually a daily alias. These peaks are identified 
for various candidate periods in Tables [H and [H for the 
combined data set and [3l and \5\ for the [Fischer et al.l 
(|2008f l data set. We also performed the same analysis 
on the combined data set of all four instruments and 
obtained consistent results. The results are also unam- 
biguous when only the Keck data are used. 
With this new period for planet e, we fit a 5-planet Ke- 

S ileria n model to the Keck and Lick data of lFischer et al.l 
2008f ). via the Leve nberg-Marqu a rdt a lgorithm im- 
plemented in IDL by iMarkwardtl (|2009D . Following 
[Fischer et al.l (|2008[ ). jitter values of 1.5 m/s and 3.0 m/s 
were adopted for Keck and Lick data, respectively, such 
that the errors became af = cr^^^t^^ ; + 4ttor,i- The re- 
sulting model fits the data much better than previous 
results, with the same number of free parameters. Com- 
pare Table [SI and Table [71 The rms is reduced from 6.45 
ms~^ to 5.91 ms~^ (10%) and the {xlV^^ is reduced 
from 1.666 to 1.411 (15%). We conclude that we have 
determined the correct period of 55 Cnc e. 

We use an epoch chosen as the weighted average of the 
observation times. The weighting was 1/cf ; this weight- 
ing minimizes the correlation between the parameters P 
and A for each planet. We have confirmed that the rms 
and {xtY^^ we achieve using a weighted epoch, as op- 
posed to using the first data point as the epoch, is iden- 
tical in the Keplerian case. 

With such a small period, we would expect planet e to 
circularize via tidal dissipation. Of course, in the pres- 
ence of perturbations of the other planets, this expec- 
tation will not be completely fulfilled. Nevertheless, we 
also repeated the fit with the eccentricity of planet e fixed 
at zero (Table [HI). 

Fitting a self-consistent Newtonian 5-planet model, 
[Fischer et afl ([20080 obtained a (xlV^^ of 2.012 
and rms of 7.712 ms~^, significantly worse than 
their best Keplerian five-planet model. We per- 
formed our own self-consistent Newtonian 5-planet 
fit using the modified Wisdorn - Holrn an symplec- 
tic integrator ('Wisdo nifc HolmanI [1991[ ) in SWIFT 
([Levison fc Dun can 1993) • Using our newly defined 
epoch, we obtain (x^)^^^ for both candidate periods of 
planet e that are statistically indistinguishable from their 
Keplerian equivalents (Table [9| and Table [TO\i . We specu- 
late that the new epoch starts the Levenberg-Marquardt 
fit closer to the global minimum and strongly recommend 
choosing the epoch as the weighted average of the obser- 
vation times, as we have done, instead of the first obser- 
vation. We have only begun to explore the dynamics of 
this system and future work adjusting the line of sight 
inclination of the system and relative inclinations of the 
planets may result in improved fits and better character- 
ization of the dynamics of this system. 
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WthMM>mmdMm^h^MdlM^^ 




Freq (1/day) 



Fig. 17. — S pectral window function of 55 Cnc for HET data set 
IjMcArthur et al. 2004). These features, convolved with a planet's 
orbi tal f requency, cause the aliases evident in the periodogram in 
Fig. [18] 




0.355 0.360 0.645 0,650 1.355 1.360 2,355 2,360 2,365 

Freq (1/doy) 

Fig. 18. — Periodo grams of 55 Cnc for pla net e only, using only 
the data from HET llMcArthur et al.ll20M) . The top row is the 
periodogram of the data themselves. The other rows show the 
periodograms of sinusoids sampled at the times of the real data 
sets as solid lines; they also repeat the periodogram of the data as 
a gray background, for comparison. Dials above the peaks show 
the phase at each peak. The second row has a sinusoid of the 
reported frequency. The third row has a sinusoid of the new fre- 
quency. In Rows 2 and 3, each peak results from the convolution 
of the sinusoidal frequency with the features in the spectral win- 
dow function in Fig. 1171 In this data set, due to noise, neither 
noiseless candidate frequency matches the data. Note the large 
phase discrepancies between the reported frequency and the data. 
Based on this data set alone, the planet's orbital period cannot be 
unambiguously determined. 
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Fig. 19. — Spectral window function of 55 Cnc for HET data set 
comb ined with ELODIE (Nacf ct al. 2004) and Lick (Marcv ot al. 
120021) . Major features of tlie spectral window function are colored: 
red (at day~^), green (yearly feature), fuschia (daily features), 
and blue (two day~^). The corresponding aliases these features 
ca use for several candidate frequencies are indicated by these colors 
in [201 
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Fig. 21. — Spectral Wi ndow Functions of 55 Cnc for combined 
Lick and Keck data set l)Fischer et al.l [20081 ). Major features of 
the spectral window function are colored: red (at day~^), green 
(yearly feature), fuschia (daily features), and blue (two day~^). 
The corresponding aliases these features cau se f or several candidate 
frequencies are indicated by these colors in l22l 
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Fig. 20. — Periodogram of 55 Cnc for planet e only. Dials above 
the peaks denote their phase. Colors correspond to the feature in 
the window function that creates the particular alias (see Fig. ll9ll . 
with red being the candidate frequency, the green sidebands yearly 
aliases, and the fuschia and blue peaks daily and two day~^ aliases 
respectively. The top row shows the data (HET+ELODIE+Lick). 
In Row 2 the solid lines show, for these time samplings, the peri- 
odogram of a sinusoid of frequency 0.3550 day~^. For Row 3, it 
is for 0.3577 day~^. For Row 4, it is for 1.3577 day~^, our now- 
favored value. The three candidate frequencies have different types 
of aliases at different locations, allowing us to break the degeneracy. 
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Freq (1/doy) 



1.560 2.555 



Fig. 22. — Periodogram of 55 Cnc for planet e only. Dials above 
the peaks denote their phase. Colors correspond to the featu re in 
the window function that creates the particular alias (see Fig. l21l l. 
with red being the candidate frequency, the green sidebands yearly 
aliases, and the fuschia and blue peaks daily and two day~^ aliases 
respectively. Row 1 shows the data (Lick+Keck). In Row 2 the 
solid lines show, for these time samplings, the periodogram of a 
sinusoid of frequency 0.3550 day~^. For Row 3, it is for 0.3577 
day"'^. For Row 4, it is for 1.3577 day~^, our now-favored value. 
The three candidate frequencies have different types of aliases at 
different locations, allowing us to break the degeneracy. 
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55 Cnc Combined Data Set: 



TABLE 2 

EXPECTATIONS FROM THE WINDOW FUNCTION. 



Candidate 
frequency, / 


0.0028-f 


0.0028+f 


1.0000-f 


Window Function Feature 
1.0027-f 1.00004-f 


1.0027+f 


2.0028-f 


2.0028+f 


0.3550 
0.3577 
1.3577 


0.3522 
0.3549 
1.3549 


0.3578 

0.3605 

1.3605 


0.6450 
0.6423 
0.3550 


0.6477 
0.6450 
0.3577 


1.3550 

1.3577 
2.3577 


1.3577 

1.3604 

2.3604 


- 


2.3604 



Note. - 
ws 2-4 
both a 


— Along the top row arc peaks in the w^indow^ function at frequencies fs (Fig;- fT9l . Each ro^v refers to a candidate frequency /; rows 1-3 in this table 
n Fig.[|o] respectively. The cells are frequency values \f ± fg\ expected for peaks in the periodogram. If the predicted alias is consistent with a peak in 
mplitude and phase, the cell is bolded. A non-emphasized cell indicates a large discrepancy in amplitude or phase. For dashed cells, no comparison ^ 


match t 
the dat 
vas done 


nits arc 


day~l- The frequency of / = 1.3577 day~l is overwhelmingly the best match to the data. 





TABLE 3 





55 


Cnc IFischer 


ET AL.I 1120081) Data S 


ET: EXPECTATIONS FROM THE WINDOW 


FUNCTION 






Candidate 
frequency, / 


0.0028-f 


0.0028-l-f 


1.0000-f 


1.0027-f 


Window Function Feature 
1.0000+f 1.0027+f 


2.0027-f 


2.0055-f 


2.0027+f 


2.0055+f 


0.3550 
0.3577 
1.3577 


0.3522 

0.3549 
1.3542 


0.3578 

0.3605 

1.3605 


0.6450 
0.6423 
0.3550 


0.6477 

0.6450 

0.3577 


1.3550 

1.3577 
2.3577 


1.3577 

1.3604 
2.3604 


0. 6450 


0.6478 


2.3577 

2.3604 


2.3605 
2.3632 


Note. — The fo 


mat is the sarr 


le aa Tabled Features in the v 


window functio 


„ are from Fig 


ET] The candidat 


frequencies i 


-1 rows 1-3 in 


this table match 


to rows 2-4 in 



TABLE 4 
55 Cnc Combined Data Set: features in the data periodogram. 



Candidate 1 
frequency, / 


0.3550 






0.3577 






Major 
0.6450 


Data Feature 
1.3577 






2.3577 




2.3604 


0.3550 
0.3577 
1.3577 




f 

f-0.0028 
f-1.0027 






f + 0.0028 

f 

f-1 .0000 






1.0000-f 
1.0027-f 
2.0027-f 




1.0027+f 

1.0000+f 

f 






2.0027+f 
1.0000+f 




2.0027+f 
1.0027+f 


Note, — The to 
match to rows 2-4 
the cell is bolded. 
at that frequency. 


P r 

A 
Ur 


3w indicates a major 
FiK.1191 respectively, 
non-emphasized cell 
its are day^l. This 


peak seer 
If, based 
indicates 
table sho 


1 in the data near the frequencies where aliases are predicted. Each row r 
on examining the plots, the frequency creates an alias that matchc 
a large discrepancy in amplitude or phase. A blanJt cell indicates 

ws that a frequency of 1.3577 day^^ is best able to account for tt 


efers to 
s that p 
that the 
e pealcs 


a candidate frequ 
3ak in the data in 
candidate freque 
m the data. 


ency; re 
both a 
ncy doe 


ws 1-3 in this table 
rrplitude and phase, 
s not cause an alias 



TABLE 5 
55 Cnc IFischer et al] II2008D Data Set: features in the data periodogram. 



Candidate 
frequency, / 



0.3550 



0.3577 



Major Data Feature 
0.6450 0.6478 



1.3577 



2.3577 



2.3604 



0.3550 
0.3577 
1.3577 



f 

f-0.0028 
f-1.0027 



f + 0.0028 

f 

f-1 .0000 



1.0000-f 
1.0027-f 
2.0027-f 



1.0027-f 



2.0054-f 



1.0027+f 

1.0000+f 

f 



2.0027+f 
1.0000+f 



2.0054+f 
2.0027+f 
1.0027+f 



; Table^ Candidate fre. 



5-8 in Fig.[2T] 1 
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TABLE 6 

55 CnC KePLERIAN radial velocity fit, Pe = 2.8 DAYS.'' 





K 


M sin i 


P 


a 


e 


LJ 


A 


Vl 


Vk 


Y^ 


N 


(x2)^/^ 


rms 




ms^"'" 


M,/„p 


days 


AU 




deg 


dog 


ms ^ 


ms ^ 








ms-i 


e 


5.2(2) 


0.0346(16) 


2.81705(5) 


0.0382(3) 


0.066(48) 


238.(41) 


86(14) 














b 


71.3(3) 


0.824(3) 


14.65164(11) 


0.1148(8) 


0.014(4) 


135.(15) 


327.4(10) 














c 


10.0(2) 


0.167(4) 


44.349(7) 


0.2402(17) 


0.09(3) 


66.(17) 


312(7) 














f 


5.3(3) 


0.148(9) 


259.7(5) 


0.780(6) 


0.40(5) 


182.(9) 


308(14) 














d 


46.9(4) 


3.84(4) 


5191.(53) 


5.76(6) 


0.015(9) 


223.(33) 


201(4) 


6.8(6) 


5.9(7) 


813.2 


27 


1.666 


6.45 



Note. — Data 
hould mi 



the 



■ et al.l l200a V T 



( the 



riic following gravitat 
a,r maBscs. Formal 
irdinatcs. as rccoir 



ck and Keck data presented by [Fit 

bion in the errors between F* and A for each planet. 

al constants were used: GMq = 0.0002959122082856, 

from the Levenberg-Marquardt algorithm are given in j 

d bv lLee h Pe^I^ HMl) . 



[JD 2453094.762), which 



tio of the 
entheses, 



Jupiter = 1047.35. Tlic : 
g to the final digit(s). Ma 



TABLE 7 
55 Cnc Keplerian radial velocity fit. 



: 0.74 DAYS.' 





K 


M sin i 


P 


a 


e 


Ld 


A 


Vl 


Vk 


X^ 


N 


ixl)'^' 


rms 




ms^-'- 


Mjup 


days 


AU 




deg 


deg 


ms ^ 


ms -"- 








ms"-"- 


e 


6.2(2) 


0.0261(10) 


0.736539(3) 


0.01564(11) 


0.17(4) 


177.(13) 


126(2) 














b 


71.4(3) 


0.826(3) 


14.65160(11) 


0.1148(8) 


0.014(4) 


146.(15) 


139.7(2) 














c 


10.2(2) 


0.171(4) 


44.342(7) 


0.2402(17) 


0.05(3) 


95.(28) 


90.(2) 














f 


5.1(3) 


0.150(8) 


259.8(5) 


0.781(6) 


0.25(6) 


180.(12) 


36(4) 














d 


46.6(4) 


3.83(4) 


5205.(54) 


5.77(6) 


0.024(10) 


192.(16) 


222.7(8) 


6.7(5) 


6.5(6) 


583.1 


27 


1.411 


5.91 



Note. — Data 
linimize the coi 



5 betwei 



cnted bv lFisch~ 
and A for each pi 



ionnooi). T 



net. For planet c 



of the 



)bscrvatioii times (JD 2453094.762), which should 
:hofTtr[,ID] = 2453094.728(10) + EX 0.736539(3). 



TABLE 8 
55 Cnc Keplerian radial velocity fit. Pe 



: 0.74 DAYS,ee 





K 


M sin i 


P 


a 


e 


Ul 


A 


Vl 


Vk 


X^ 


N 


ixiy^' 


rms 




ms~^ 


Mj„p 


days 


AU 




deg 


deg 


ms ^ 


ms ^ 








ms""*" 


e 


6.1(2) 


0.0258(10) 


0.736540(3) 


0.01564(11) 


0.000(0) 


0.(0) 


126(2) 














b 


71.4(3) 


0.825(3) 


14.65158(11) 


0.1148(8) 


0.012(4) 


147.(17) 


139.7(2) 














c 


10.3(2) 


0.172(4) 


44.341(7) 


0.2402(17) 


0.06(3) 


99.(23) 


90.5(15) 














f 


5.0(3) 


0.150(8) 


260.0(5) 


0.781(6) 


0.13(6) 


180.(21) 


37(3) 














d 


46.7(4) 


3.83(4) 


5214.(54) 


5.77(6) 


0.029(10) 


189.(14) 


222.6(8) 


6.8(5) 


6.3(6) 


598.1 


27 


1.429 


5.89 



Note. — Data 
planet the orbit i 



ck data pn 
s held at zc 



ited by [FT 
For plane 



m nnna ). 



-,^^[JD] ^ 2453094.762 Because tidal dissipation has most likely nearly 
edict a transit epoch of T^j.[JU] = 2453094.688(4) + E X 0.736540(3). 
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4. DISCUSSION 

4.1. Summary of Approach 

Aliases result from a convolution between a true phys- 
ical frequency and the spectral window function, which 
is created by gaps in the data set due to observational 
constraints. Our method harnesses features in the win- 
dow function to distinguish aliases from true frequencies. 
For a given frequency / and window function peak /s, 
aliases will occur at |/ ± /s|, where /^ is a feature in the 
window function. In the ranges where we expect major 
aliases to occur, we compare the phase and amplitude 
of aliases predicted by a sinusoid of the candidate fre- 
quency sampled to the data, with other known planets 
subtracted off beforehand. We judge whether the "pat- 
tern" of the predicted aliases matches the data: for exam- 
ple, yearly aliases appear as sidebands of the candidate 
frequency while daily aliases often appear as a doublet 
caused by the sidereal and solar day. If all the aliases 
match in amplitude, phase, and pattern, we can be con- 
fident that we have found the true orbital period. If 
there are discrepancies and the aliases of none of the can- 
didate frequencies match the data, we know that noise 
prevents us from definitively determining the true period 
and that follow-up observations are necessary. Misunder- 
standings about aliases have previously led to incorrect 
identification of planet's orbital periods, a key param- 
eter in defining the planet's properties, as well as the 
dynamical behavior of the planets in the system. We 
have corrected common misconceptions, including that 
aliases always appear near the frequency of peaks in the 
window function, that any frequency above 1 cycle/day 
is necessarily an alias, and that aliases will appear if the 
data are scrambled or if the true frequency is subtracted 
out. 

4.2. Summary of Results 

For two systems, we confirmed previous distinctions 
between alias and true frequency. The period of G J 876 d 
is indeed 1.94 days, not 2.05 days. The period of 
HD 75898 b is indeed 400 days and the periodogram peak 
at 200 days is indeed an alias, not a second planet or 
eccentrici ty harmonic, the alter native explanations pro- 
posed bv lRobinson et all (|2007[ ). 

For two other systems, we determined that the data are 
too noisy to allow us to definitely distinguish between 
alias or true frequency. According to our analysis, it 
remains unclear whether the period of Gl 581 d is 67 
days or 83 days; even a period of 1 day cannot be ruled 
out. It also remains unclear whether HD 73526 contains 
two planets with orbital periods 187.5 and 376.9 days, 
locked in a 2:1 resonance, or whether one of the periods is 
actually 127 days. Further observations of these systems 
are required, preferably at times that reduce the aliasing. 

For a final pair of systems, we determined the reported 
orbital period was incorrect, due to mistaking a daily 
alias for the true frequency. According to our analysis, 
the orbital period of HD 156668 b is actually 1.2699 days, 
not 4.6455 days. The orbital period 55 Cue e is 0.7365 
days, not 2.817 days. The standard, general-purpose 
software Si qSpec mentioned in the introduction (Reegen 
l2007l l2010t ) agrees with our orbital period distinctions 
(we used the parameters: depth=2, par — 0.2 and par = 
0.5, and a frequency upper limit of 2 day— 1). 



4.3. Implications for 55 Cnc e 

What are the implications of an updated period for the 
innermost planet of 55 Cnc? 

First, it dramatically lowers the effective noise when 
determining t he pa rameters of the planetary system. 
iFischer et al.l (|2008[ ) reported independent Keplerian fits 
with rms of 6.74 m/s, and a self-consistent dynamical fit 
with rms of 7.712 m/s. Our Keplerian fit achieves rms of 
5.91 m/s, and our self-consistent coplanar dynamical fit 
achieves rms of 5.96 m/s. By adjusting the inclination of 
the system relative to our line-of-sight and the planets' 
mutual inclinations, an even better self-consistent might 
be possible. Therefore perturbations might be directly 
detected via a lower rms when interactions among the 
planets are included, and the architecture of the system 
further constrained. We have just begun exploring this 
avenue. 

Second, 55 Cnc e itself can now be searched for tran- 
sits at the new period, with high a priori probability 
of ^ 25%. Given the period and phase of the radial- 
velocity signal, we report predicted transit epochs in Ta- 
bles [7] and [H The predictions differ because the lat- 
ter assumes zero eccentricity, and the formally signifi- 
cant value of Cc matters. Nevertheless, folding the sys- 
tematic uncertainty related to eccentricity into the pre- 
dicted transit time, we still can predict transit times good 
to (Tt ~ 1 hour in 2010. This search can be accom- 
plishe d simply by folding the photometric data reported 
by Fi scher et al.l ([2008) at the new ephemeris. Gregory 
Henry (priv. comm.) has made such a search, finds no 
positive signal, and constrains putative transits in the 
period range 0.7 — 0.8 days to a depth < 0.7 mmag, or 
> 2.6 i?©. Earth-composition models of super-Earths 
predict a radius ^1.9 i?® (Valencia et al. 2006), so a 
search at higher precision is certainly worthwhile. 

Third, even apart from a transit, this super-Earth must 
be very h ot, as it is very cl ose-in to a solar-type star. 
Following iLeger et al.l (|2009f ) , we find that the substel- 
lar point could be up to 2750 K, if the insolation is ab- 
sorbed then reradiated locally. We would naively expect 
that the enormous radiation this planet takes in woul d 
evaporate any atmosphere (e.g., iJackson et al.l l2010( ). 
Moreover, the host star is also very bright as seen from 
Earth. Therefore it might be useful to look for its phase 
curve with Spitzer , to d etect or rule out an atmosphere 
(jSeager fc Demingll2009f ). Another attractive possibility 
is probing a ma ,gma ocean, which ruay exist because o f 
the irradiation (jGelman et al.ll2009HGaidos et al.ll2010D . 
but this may require transit measurements. 

Fourth, the presence of the other 4 planets surely in- 
jects a non-zero eccentricity into this tidally-dissipating 
planet. Its expected value remains to be calculated, but 
will likely be on the order of 10~*. This forced eccen- 
tricity could stimulate co nsiderable geologic activity — 
it might be a "super-Io" (jBarnes et al.ll2010r) . 

4.4. Observational Strategies for Mitigating Aliases 

Can aliases be prevented or mitigated by the choice of 
observation times? Constraints on when the star is vis- 
ible at night necessarily result in gaps in the data that 
cause aliases. However, we encourage ob servers to en- 
gage in "window carpentry" ()Scarglelll982[ ) by observing 
the star during the greatest span of the sidereal and solar 
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day possible, not just when the star transits the meridian. 
Unfortunately, observing stars as they rise and set poses 
a challenge for observers, who minimize slew time^ and 
thus maximize the number of stars observed per night by 
observing at the meridian for the majority of the night. 
This observing strategy (Fig. [SJd) results in strong daily 
and yearly aliases (ex. Fig. [T31) . Another strategy is to 
start in the west and gradually move east over the course 
of the night (Fig. [S]:), observing as much of the sky as 
possible. This strategy reduces yearly aliases but side- 
real daily aliases remain strong (ex. Fig. [TS] and I^Tj) . 
To reduce sidereal daily aliases, we recommend the fol- 
lowing procedure. Start the telescope somewhere west 
of the meridian (randomized from night to night) and 
move east to cover half the sky over the course of half 
the night. Then make one large slew to the place the 
telescope started and re-observe the same portion of the 
sky. Some stars will gain the advantage of being ob- 
served twice in one night. Moreover, when the data are 
folded at the mean sampling period, they still show some 
variety in phase of observation, which is needed to re- 
duce window function peaks and de-alias candidate pe- 
riods. However, another consideration is that at higher 
air mass, both the extinction is greater and the seeing is 
worse. The increased atmospheric attenuation means a 
longer integration time is required, reducing the number 
of stars that can be observed, while the seeing increases 
the measurement errors. For a particular set of stars, ob- 
servers can work out a slew pattern that will maximize 
the number of stars observed while minimizing aliasing. 
[Saunders ct al. (200§) present a clever method for de- 
termining the optimal sampling when period searching 
using satellite telescopes or a longitude-distributed net- 

^ If tho slew time exceeds the read-out time, fewer observations 
may be gathered per night. However, the wise spacing of observa- 
tion times can more than make up for this through disambiguation 
of ahas frequencies using fewer data points. 



work that can observe continuously. Unfortunately this 
strategy is impractica l to implem ent using a single tele- 
scope on the ground. iFordI (|2008[ ) presents useful adap- 
tive scheduling algorithms for observing multiple targets 
that can be parametrized to reduce aliasing. 

We suggest taking advantage of any unusual time win- 
dows: for example, the rare granted dark time or time 
at the beginning or end of another observer's night. Ob- 
servers focusing on a large group of stars can determine 
which star would most benefit from this unusual time 
by calculating the window function with the new obser- 
vation times added or, in the case of a planet with two 
candidate periods, determining for which system the ob- 
servation times would best distinguish between two can- 
didate orbits. We also suggest that it would be beneficial 
to observe stars using telescopes in two or more locations 
at different latitude and longitudes (ex. Fig. [T9| . 

At the stage of data analysis, we encourage the use of 
our method to distinguish true frequencies from aliases, 
crucial for the correct characterization of the planet. As 
astronomers push to observing lower mass planets and 
modeling planets near the noise limit, they cannot as- 
sume that the highest peak in the periodogram - or even 
the best Keplerian fit - corresponds to the true orbital 
period. Only by harnessing features in the window func- 
tion to compare the amplitude, phase, and pattern of an 
assortment of predicted aliases to the data can we distin- 
guish the planet's true orbital frequency - or determine 
that more observations are needed. 

4.5. Conclusion 

Knowing a planet's correct orbital period is essential 
for accurately characterizing it. By Kepler's law, the 
planet's distance from the star increases as its orbital pe- 
riod increases. Therefore the planet's orbital period sets 
its temperature: too hot, too cold, or just right for life. 
The planet's inferred mass, as calculated from the radial 
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velocity amplitude, increases as the period decreases - a 
closer planet needs less mass to exert a given force on 
the star - so a difference in orbital period may be the 
difference between an Earth analog and a super-Earth. 
In the case of multi-planet systems, the spacing of the 
planets determines their mutual interactions: therefore 
a difference in orbital period may be the difference be- 
tween a precariously placed planet and one locked deep 
in a stabilizing resonance. The signal of a planet's ec- 
centricity is contained in the harmonics of the planet's 
orbital period: therefore a difference in orbital period 
may be the difference between a planet that formed in 
situ and a planet violently scattered, a calm planet that 
has long been tidally circularized or a planet erupting 
with volcanoes due to tidal dissipation. But periods that 
correspond to totally different worlds are only subtly dis- 
tinguishable in the radial velocity signal. Such are the 
machinations of aliases. 

Through our method, astronomers can confirm a 
planet's orbital period or determine that noise prevents a 
definitive distinction. In the latter case, follow-up obser- 



vations taken according to the suggestions above should 
eventually allow the true period to be determined. Iron- 
ically, Earth's own rotational and orbital period make it 
challenging to uncover the orbital period of other worlds, 
particularly Earth analogs. But by better understanding 
of digital signal processing, we can mitigate the deleteri- 
ous effects of the inevitable sunrise and starset. 
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